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Abstract 

Partitioning and transport of water and small solutes into and through nanopores is important to a variety 
of chemical and biological processes and applications. Here we study water structure in negatively charged 
model cylindrical [carbon nanotube (CNT)-Uke] nanopores as well as the partitioning of positive ions of 
increasing size {Na'^, K^, and Cs^) into the pore interior using extensive molecular dynamics simula- 
tions. Despite the simphcity of the simulation system - containing a short CNT-hke nanopore in water 
carrying a uniformly distributed charge of Qpcn-e = —ne surrounded by n (= 0, . . . , 8) cations, making the 
overall system charge-neutral - the results provide new and useful insights on both the pore hydration and 
ion partitioning. For n = 0, that is, for a neutral nanopore, water molecules partition into the pore and form 
single-file hydrogen-bonded wire spanning the pore length. With increasing n, water molecules enter the 
pore from both ends with preferred orientations, resulting in a mutual repulsion between oriented waters at 
the pore center, and creating a cavity-like low density region at the center. For low negative charge densities 
on the pore, the driving force for partitioning of positive ions into the pore is weak, and no partitioning is 
observed. Increasing the pore charge gradually leads to partitioning of positive ions into the pore. Inter- 
estingly, over a range of intermediate negative charge densities, nanopores display both thermodynamic as 
well as kinetic selectivity toward partitioning of the larger and Cs~^ ions into their interior over the 
smaller Na'^ ions. Specifically, the driving force is in the order > Cs'^ > Na'^, and K'^ and Cs"*" 
ions enter the pore much more rapidly than iVa+ ions. At higher charge densities, the driving force for 
partitioning increases for all cations - it is highest for if + ions, and becomes similar for Na'^ and Cs'^ 
ions. The variation of thermodynamic driving force and the average partitioning time with the pore charge 
density together suggest the presence of free energy barriers in the partitioning process. We discuss the 
role of ion hydration in the bulk and in the pore interior as well as of the pore hydration in determining the 
barrier heights for ion partitioning and the observed thermodynamic and kinetic selectivity. 
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I. INTRODUCTION 



Understanding thermodynamics and kinetics of partitioning of small molecules into nanoscopic 
pores or confined spaces of molecular dimensions is important to a wide range of chemical and 
biological processes and applications llli|2i|3i|4|]. The pores of interest can be smooth cylindrical 
(carbon nanotube-like) Jsl] or slit-pores (between flat surfaces) or may have complex geome- 
tries as observed in biological pores 171,180, in activated carbon, and in other nanoporous materials 



. When exposed to solution, such materials can selectively extract and/or allow transport of spe- 
cific solute molecules into and through their pores, depending on the pore-solvent, pore-solute, and 
solute- solvent interactions. Successful bottom-up design of such selective systems can be aided 
by fundamental understanding of the interplay of these interactions in idealized model systems. 
Molecular dynamics simulations are an excellent tool to build and simulate such model systems 
aimed at fumdamental understanding. 

Here we are interested in partitioning of simple spherically symmetric ionic solutes into interior 
of charged hollow cylindrical pores of well defined dimensions in the nanometer range. Such a 
model system is relevant not only to electrostatically driven molecular separations or nanofluidic 
applications, but also to development of energy storage devices, such as carbon supercapacitors, 
that use ion adsorption on the surface of highly porous materials to store charge. A recent exper- 
imental study by Chmiola et al. [^1^ shows that carbide-derived carbon materials with pore sizes 
from 0.6 to 2.25 nm when used as a negative electrode, can partition tetraethylammonium cations 
(~0.68 nm in diameter) into the pore. The closer approach of partially dehydrated cations to the 
nanoporous electrode surface leads to significantly enhanced capacitance of these materials. 

Our simulation system contains a nanopore (modeled by a carbon nanotube-like cylindri- 
cal pore) in water carrying a uniformly distributed charge of g^ore = —ne surrounded by n 
(= 0, . . . , 8) cations (either Na'^, K^, or Cs+), making the overall system charge-neutral. We 
perform extensive nonequilibrium simulations focused on the kinetics of partitioning, which show 
that negatively charged nanopores can be selective toward partitioning of the larger and Cs"*" 
ions and can exclude the smaller Na'^ ions over a range of nanopore charge densities. We com- 
plement those simulations with equilibrium free energy calculations of ion hydration in bulk water 
and in the hydrated nanopore interior. These calculations along with analysis of ion hydration 
shell fluctuations highlight the role of free energy barriers as well as equilibrium free energy of 
partitioning which together are expected to govern the rate of ion transport through selective pores. 
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II. SIMULATION DETAILS 



Simulations of ion partitioning into nanopores: We used a piece of (5,5) armchair carbon nan- 
otube (CNT)-like pore comprising 100 carbons [with the pore diameter {i.e., C-C distance) of 6.7 
A and length 1 1 A] as a model cylindrical pore. A total charge of q^ore = (where e is the 
magnitude of electronic charge) was distributed uniformly on the pore atoms, such that each atom 
carries a charge of -ne/100. The nanopore was placed in a solution containing 500 explicit water 
molecules and n cations, thus, making the overall system electrically neutral. Note that our goal 
is not to model specifically a CNT, but use its cylindrical geometry as a model for a cylindrical 
pore. A realistic model of CNT would need to consider the chemistry of rim atoms, non-uniform 
distribution of charges (in case of charged nanotubes) lllin . etc. 

The present systems contain n cations, but no explicit anions (the pore with its —ne charge 
acts as a large anion). In realistic systems, anions (e.g., chloride ions) will be present in the salt 
solution. For n = 5, we simulated the present system with excess salt {i.e., with additional few 
pairs of NaCl, KCl, or CsCl) ions. The results on partitioning of cations in this system were 
similar to the ones reported here, suggesting that the trends on partitioning reported here will be 
unaffected. 



Simulations were performed using AMBER6.0 [|12|1 . using TIP3P [|13|] model of water, and 
Lennard Jonesium (LJ) description of carbon (atom type CA) [12, Im]. Ions were represented as 
LJ spheres with ion charge placed at the center 111511 : ciNa = 2.530 h, e^a = 0.06184 kJ/mol, 
aK = 5.874 A, eK = 0.0000568 kJ/mol, and cxcs = 6.0492 A, e^, = 0.0003372 kJ/mol. Ion- 
water and ion-carbon LJ interactions were calculated using Lorentz-Berthelot mixing rules [llql . 
Periodic boundary conditions were applied and the particle mesh Ewald method [ | 17|1 was used to 
calculate the electrostatic interactions with a grid spacing of 1 A. Temperature and pressure were 
maintained at 300 K and 1 atm, respectively, using the Berendsen method UM- A time step of 1 
fs was used in all the simulations. 

Selective partitioning of cations into nanopores: A given system contains a carbon nanotube- 
like pore carrying a total charge of Qpore = —ne, water, and n cations. We monitored the number of 
cations partitioned into the nanopore as a function of time for different values of n (see Figure [l]). 
Separate simulations were performed for qpore = — le through — 8e in steps of le, for each cation. 
At the beginning of each production run (marked t = 0), each cation in the system was at least 9 A 
away from the atoms of the model nanopore. To obtain an estimate of the partitioning time (or the 
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first passage time for partitioning), each simulation was run for a sufficiently long time. For lower 
charge densities on the nanopore, the electrostatic driving force is low, and extended simulations 
of neutral and near neutral {qpore = — le ) nanopore s show that ions are excluded from the pore 



interior, consistent with previous studies [|19LI2(a |21|1. For higher charge densities on the pore, not 



only is the electrostatic driving force higher, but the ion concentration outside the tube is slightly 
higher as well (by the way the system is constructed), which likely further increases the driving 
force. Thus, with increasing n we expect ions to partition into the pore, with the time required for 
partitioning decreasing with the pore charge density. Correspondingly, for Ignore I > 2e simulations 
were performed for 24/n nanoseconds each. 



III. RESULTS AND DISCUSSION 

Kinetics of ion partitioning: Figure \T\ shows the number of cations partitioned into the pore 
for various charge states of the pore observed in different nonequilibrium simulation runs. For 
Qpore = — 2e, the charge density on the nanopore is sufficiently small and none of the ions partitions 
into the pore over the timescale of 12 ns. It is only for qpore = — 5e we observe partitioning of one 
Na~^ ion into the pore in the specific set of runs shown in Figure [U In contrast, Cs^ and ions 
partition into the pore for Ignore I > 3e. At higher values of \qpore\, we observe partitioning of the 
second and third (in case of K^) cation into the charged nanopore. To obtain more quantitative 
estimates of the kinetics of partitioning, we performed 10 independent nonequilibrium simulations 
for each ion type starting with different initial configuration for Qpore = — 4e, — 5e, and — 6e. 
Each of these simulations were run until the first cation partitions well into the pore, which for 
Qpore = — 4e was longer than 25-30 ns each for several runs with Na'^ ions. We define ti as 
the waiting time for the first cation to partition into the nanopore. Table IJ lists average, (ri), and 
standard deviation, a^, of the waiting time distribution obtained from 10 simulations. Although 
a larger number of simulations may be needed to obtain more accurate estimates of waiting time 
distributions, trends in our data are already clear and interesting. For example, for qpore = — 4e, 
Na~^ ion takes ~25 times longer to partition into the pore (~18 ns) compared to that for and 
Cs~^ ions (~0.7-0.8 ns). As the charge on the nanopore is increased, partitioning of Na'^ occurs 
over shorter timescales, and correspondingly, the ti values for both and Cs+ ions are small 
and decrease somewhat further. Thus, from kinetic perspective, over a range of charge densities, 
the nanopore is selective toward partitioning of the larger ions and Cs^ over the smaller ion 
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Na^, and increasing the driving force significantly (i.e., increasing |gpore|) appears to decrease 
that selectivity somewhat. To understand the ion partitioning at a more fundamental level, we 
investigate below the structural and thermodynamic aspects of the partitioning process. 



il 



ing and emptying transitions of water into 



2M- For the carbon- water LJ parameters used 



Water structure in neutral and charged pores: \ 
neutral nanotubes have been studied previously 121 
here, water molecules partition into the tube and form a single-file hydrogen bonded wire, in which 
each water molecule donates and accepts one hydrogen bond to and from its neighbors on the left 
and right (Figure [2k)- The length of the tube is such that only four water molecules completely fit 
inside the tube. However, the hydrogen-bonding interactions of the boundary water molecules with 
those inside the tube are sufficiently strong. As a result, one water molecule on each side of the 
tube also maintains its position and orientation consistent with those in the tube. Correspondingly, 
the water density profile along the tube axis shows six clearly defined peaks (Figure [2b)- 

Water molecules in the hydrogen-bonded wire have one of their OH bond vectors aligned ap- 
proximately with the pore axis and pointing in the positive Z direction. Properties of such hydro- 
gen bonded wires have been studied especially for their ability to transport protons at fast rates 
[[2311 ■ The dipole vector OM of each water molecule makes an angle of ~35 degrees with the pore 
axis \_{cos9) ^ 0.82] (Figure [Ifc). We note that the mirror image of the chain with all OH vectors 
pointing in the negative Z direction is equally likely due to two-fold symmetry ^2!^ . 

The effects of charging the nanopore on the structure of water inside the pore are interesting. 
For Qpore = — 2e, that is, when the charge on each carbon is — 0.02e, water molecules enter the 
pore with both of their hydrogen atoms pointing inward and the HOH planes roughly parallel to the 
pore axis (Figure [2h) L22]. As the chains from the two ends meet at the center, there is significant 
electrostatic repulsion between them. As a result, water density is depressed at the center of the 
pore as indicated by the significant reduction in the heights of the two central peaks (Figure [2b). 
Water is effectively expelled from that region, creating a molecular scale void or a cavity. This 
expulsion in similar to that observed recently in simulations of charged plates in water |[24ll . The 
other peaks are also less well defined. As the charge on the nanopore is increased, the stronger 
water-nanopore interactions lead to a partial filling of the internal cavity. In the interior of these 
charged nanopores, water molecules flip the orientation as we traverse the pore from left to right 
as indicated by the {cosO) profile (Figure[2fc)- A similar change in orientation of water molecules 
is observed in the aquaporin protein channel, vvhich has positively charged regions at the opening 



and the exit to block transport of protons [[7, 
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The low density region at the center of a charged pore is only a few angstroms wide and presents 
an ideal location for cations partitioned into the pore. The lower water density in the cation sized 
region ensures lower repulsive interactions without the loss of attractive interactions with vicinal 
waters. Indeed, analysis of density distribution of cations confirms that the most favorable po- 
sition for the partitioned cations is at the pore center. The cation-water electrostatic interactions 
are strong and affect orientations of the vicinal waters. Figure |3] shows the effects of competi- 
tive water- water, water-nanopore, and water-cation interactions on water orientations. When the 
charge on the cation at the pore center is turned off, e.g., for K'^Fore^^^ case, water orientations 
are dictated solely by the water-nanopore electrostatic interactions. However, when the charge 
on the potassium ion is turned on, i.e., for the K+^Pore^^*^ case, the water-potassium interactions 
dominate and flip the orientations of two vicinal water molecules such that their oxygens point to- 
ward the ion. Interestingly, the repulsion between these two water molecules in ion's shell and 
others located at the pore ends leads to a low density area at the pore entry (see Figure [3t). When 
the nanopore charge is increased to — 5e, the water-nanopore interactions become comparable to 
water-potassium interactions, and only the two water molecules vicinal to the ion point their 
oxygens toward the ion, and the pore entry region gets filled again (Figure [3]1). 

Thermodynamics of ion partitioning: The free energy of ion transfer from bulk water to 
nanopore interior quantifies the thermodynamic driving force for ion partitioning. The excess 
free energy of hydration, /i^^^, can be divided into two parts: fifj^^ = ji^^j + /x^fg, where is 
the excess chemical potential of hydration of electrically neutral LJ solutes, and /i^fg is the free 
energy of charging that solute to its final charge state (+le here). We used test particle insertions 



2711 of neutral LJ solutes in bulk water and at the center of a 



charged hydrated nanopore to 



obtain the contribution. As described previously in detail [|28l.l29(]. calculation of the /Xgfg was 
done using the cumulant expansion method which requires two simulations, one in the uncharged 
and another in the fully charged state of the solute in bulk water as well as in the nanopore in- 
terior. Such calculations were performed for all three cations for nanopore charge of — 2e, —he, 
and — 8e. For reference, we also calculated /i^j and /i^fg contributions in a system containing a 
charged nanopore in vacuum, in the absence of water molecules. Comparison of these values to 
those in the hydrated nanopore systems provides insights into the role of water in ion solvation, 
especially in the confining pore interior region. 

Table In] lists values and /i^fg for Na^ , , and Cs+ ions in bulk water, in the hydrated 
nanopore interior, and at the center of nanopore in vacuum for different charge states of the 
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nanopore. In bulk water, the is smallest (i.e., most favorable) for A^a+ ion, and increases 
with the ion size. This is expected because the yu|j is dominated by repulsive interactions (or 
the cavity formation process). Compared to that in bulk water, value in the interior of the 
hydrated nanopore with — 2e charge is smaller for Na'^ ion (expected), somewhat higher for 
ion, and almost doubles of the Cs~^ ion (unexpected). This is surprising because we expect that 
the presence of a low density region in the hydrated nanopore will reduce, and not increase the 
/ifjj value for all ions. The trends become clear, however, when we consider fxfj calculated for 
ions in a nanopore in vacuum. They indicate that the interior pore of the (5,5) nanopore used here 
is small; it is large enough to accommodate Na'^ ion comfortably, but overlaps slightly with the 
ion, and significantly with the Cs+ ion. 
The advantage arising from the presence of a cavity at the center of the hydrated nanopores 
becomes clear when we compare the sum of in bulk water and in bare nanopore to that in 
hydrated nanopore. For all ions, that value in hydrated nanopores is smaller than the sum of 
bulk water and bare nanopore values. Only for the highest charge density on the nanopore, the 
two approach each other as the central cavity gets gradually filled with water. Thus, based on 
differences between /i|fj in bulk water and in the hydrated nanopore, we would expect the hydrated 
nanopore interior to be most favorable for the smallest Na~^ ion and least favorable for the largest 
Cs+ ion. This is contrary to the kinetic behavior observed in Figure [T] To understand the role of 
electrostatic interactions and of hydration, below we focus on differences in nlf^ contribution for 
the three ions. 

Table In] shows that in bulk water phase, fj^lf^ is large and negative (i.e., favorable) for all ions. 
It is most favorable for Na'^ ion, equal to -420 kJ/mol, and decreases in magnitude as the ion 
size increases, consistent with previous calculations [|3ol] . In vacuum, the electrostatic interaction 
between cations and the negatively charged nanopore is highly favorable and is almost identical 
for all ions. Minor differences in the energy values in vacuum arise due to slight differences in the 
location of ions. Ions are placed at their most probable locations obtained from ion partitioning 
data in the hydrated charged nanopores. The charging free energy of ions in the hydrated charged 
nanopores is also large and negative. However, its magnitude is smaller than that in vacuum 
calculations. In the charged and hydrated nanopore, water molecules from either side point their 
hydrogen atoms toward the center. Charging of a neutral LJ solute located at the center to +le 
charge partially reorients those water molecules to more favorable ion-water configurations. This 
effect can be viewed as "competitive solvation", where the cation and the pore both compete for 
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the same water molecules for their solvation. As a result, although the net value of /i^fg for ion 
solvation is negative, it is smaller in magnitude compared to that inside the nanopore in vacuum. 
Sodium ion, with its small size, and higher charge density is able to fit in, as well as orient vicinal 
water molecules in the pore such that the overall reduction (compared to that in vacuum) is smaller 
compared to that for potassium and cesium ions. 

The sum, fxf^^ = + fxlf^, determines the overall driving force for ion partitioning into the 
pore. The difference of ion hydration free energy in the hydrated nanopore interior and in bulk 
water is the water-to-nanopore transfer free energy, and is listed in Table UlI] as well as shown in 
Figure |4] for different charge states of the nanopore. For qpore = — 2e, the transfer free energy for 
all ions is positive (unfavorable) consistent with the observation from Figure[T]that none of the ions 
partitions into the pore over 12 ns timescale. For Na~^ and ions, the electrostatic contribution 
is unfavorable; that is, the overall ion-water-nanopore interactions do not compensate for the loss 
of ion-water interactions in the bulk water. In contrast, for Cs+ ion, it is the Afif^j contribution 
that is unfavorable as expected from the somewhat larger size of that ion relative to the nanopore 
internal diameter. 

Increasing the charge density on the nanopore increases the strength of cation-nanopore inter- 
actions, making the pore interior increasingly favorable for all three cations. At a sufficiently high 
value of nanopore charge (e.g., Qpore = — 5e), the ion-nanopore electrostatic interactions are large 
and negative for all ions indicating that it is thermodynamically favorable for all three types of 
ions to partition into the pore. The free energy varies approximately linearly with the nanopore 
charge, and the threshold value of charge at which the transfer free energy changes from positive 
to negative can be estimated from Figure IH That threshold charge value is near — 2e for and 
Cs+ ions, and somewhat larger (~ — 3e) for the Na~^ ion. That is for nanopore charge of slightly 
larger than — 2e, partitioning of and Cs+ ions will be quite favorable, whereas that of Na~^ 
ion unfavorable, providing a significant thermodynamic selectivity for partitioning of those larger 
ions over the smaller A^a+ ion. At high nanopore charge densities, although partitioning of all ions 
into the nanopore is favorable, the ions appears to benefit from its optimal size (compared to 
the nanopore diameter) and corresponding hydration free energy in water. The balance of LJ and 
electrostatic interactions for ion in bulk water and in the nanopore interior are such that its 
partitioning is most favored at higher nanopore charge densities. 

Based on the large thermodynamic driving force at charge densities greater than Qpore = — 3e, 
we would expect partitioning of all three ions into the charged nanopores. Corresponding average 
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partitioning times are listed in Table HI It is clear that the timescale for partitioning of A^a+ ions 
for qpore = — 4e is significantly larger than that for and Cs+ ions, which partition over sub- 
nanosecond timescales. The expected partitioning of Na+ in pores with total charge of — 3e would 
take even longer time, well beyond the scope of present simulations. Thus, the thermodynamic 
and kinetic data suggest the existence of a free energy barrier for the ion partitioning process, the 
height of which is different for different ions. Below we discuss possible physical origins of that 
barrier. 

Possible origins of the free energy barrier based on the dehydration of ions: The diameter of 
the nanopore is small, just sufficient to accommodate the cations. As seen in Figure [3l only two 
water molecules are available for direct hydration of ions partitioned into the pore. In contrast, the 
number of water molecules in the hydration shell of cations in bulk water is larger, approximately 
5.3, 7, and 8, for Na^ , K and Cs+ ions, respectively, consistent with previous simulation and 



experimental studies 113 lLl32l l33l. 13411. Ions thus undergo significant dehydration as they partition 



into the nanopore. As the ion approaches the end of the nanopore from outside, the ion-nanopore 
interactions partially compensate for the loss of ion-water interactions. One would require nu- 
merous simulations of the_|)artitioning process or apply a combination of methods such as the 
transition path sampling [35] and umbrella sampling [|36l. l37|l to obtain detailed structural insights 
into the transition state, the appropriate reaction path, and the free energy along that path. We have 
not performed that analysis here. However, monitoring fluctuations of the number of hydration 
shell water molecules for cations in bulk water provides qualitative insights into the origin of the 
free energy barrier. 

Figure [5] shows the potential of mean force, -h\[p{Nhy(i)], obtained from probability, p{Nhyd), 
of observing Nhyd, number of water molecules in the hydration shell of three different cations. 
Nhyd fluctuates between 4 and 7 for Na'^ ion, with an average value of 5.3. The hydration shells 
of and Cs+ display higher fluctuations, with Nhyd ranging between 4 to 10 for and 4 
to 12 for Cs+, respectively. The trends in these fluctuations are consistent with ion-water radial 
distribution functions shown in Figure [51 which show tighter binding of water molecules to Na'^ 
compared to that for the larger or Cs^ ions. The hydration number of 2 as required in the 
nanopore interior is never observed spontaneously in simulations in bulk water, as the free energy 
of such dehydration is rather high. However, the estimated value of free energy for N^yd = 2 in 
bulk water, approximately equal to 10A;T for Cs+, llkT for K^, and over 25A;T for Na'^ ion 
is consistent with the significantly slower partitioning of Na'^ compared to and Cs^ ions. 
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Despite this consistency, it is important to note that the dehydration free energies do not equal the 
free energy barrier for partitioning of these ions. The favorable ion-nanopore interactions assist 
the dehydration of cations, therefore reducing the barrier height. This is supported by the fact 
that the kinetic selectivity decreases with the increasing of the nanopore charge density. It will be 
interesting to investigate the exact nature of the barrier and especially the role of water structure at 
the pore opening in future simulation studies. 



IV. CONCLUSIONS 



Nanoporous materials with well defined pore sizes can provide excellent means to extract or 



39L 



41 



421]. However, solvent 



separate small solutes based on steric or shape effect alone [|38L 
can play an important role through differential solvation of solutes in the bulk and in the pore 
interior thereby providing selective partitioning of specific solutes of interest from a mixture of 
similarly-sized solutes. We studied the partitioning of cations Na~^, and Cs+ ions from bulk 
water into negatively charged model cylindrical pores with increasing pore charge density using 
MD simulations. Neutral and near neutral nanopores exclude ions from their interior. However, 
over a certain range of negative charge density on the pore, the nanopores display selectivity toward 
partitioning of the larger cations and Cs+ over that for the smaller Na'^ ion. Nonequilibrium 
kinetic simulations show that the partitioning is significantly slower for Na'^ compared to that 
for and Cs+ ions, especially for lower charge densities on the pore. Thermodynamic and 
kinetic observations collectively suggest the presence of a barrier for partitioning of cations into 
the nanopore interior. Such a barrier is expected to slow down the flow rates of ions through 
nanopores. We found that trends in free energy of dehydration of ions in bulk water are consistent 
with the partitioning kinetics. 

Although our studies here are focused on a CNT-like model pore, the results are relevant to stud- 
ies of carbon nanotubes focused on solute partitioning and transport 



and to complementary theoretical analyzes 147 . 
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49, 



50,51 



52L 



43 



45 



531] . Previous studies of neutral 



hydrated carbon nanotubes show that water flows through their greasy interior in a frictionless 
manner and the flow rate is limited by the entry and exit barriers [|l9ll . The low density region 
at the center of charged nanopores observed here would present an internal barrier for water 
transport. It will be interesting to explore to what extent such internal barriers will affect the flow 
rates of water through nanopores. More broadly, manipulating the balance of various interactions 
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5611 provides a mechanism for 



through design of nanosystems and functionalization [|54l . 155 
selective partitioning, separation, or transport of solutes of interest from a mixture of molecules 
of similar size and shape. Lastly, the partitioning of cations into negatively charged pores of 
subnanometer dimensions is qualitatively consistent with the picture emerging from experimental 
studies of carbon supercapacitors jioll . Simulation studies of more realistic systems along these 
lines may aid better design of high energy storage devices. 
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TABLE I: The average, (ri), and standard deviation a-j-^ of the first passage time (in ns units) for partitioning 
of iVa+, K'^, and Cs'^ ions into negatively charged nanopores at varying charge densities. 





nanopore ( 


-4e) 


nanopore ( 


-5e) 


nanopore (— 6e) 




(n) 




(ri) 




(n) o-ri 




17.46 


11.79 


3.09 


2.06 


0.92 0.65 


K+ 


0.66 


0.50 


0.72 


0.52 


0.14 0.04 




0.78 


0.41 


0.36 


0.18 


0.34 0.34 



TABLE 11: Contributions from LJ (/xfjj) and electrostatic interactions (/Xgfg) to the free energy of hydration 

of cations in bulk water, in hydrated nanopores, and in nanopores in vacuum (the last three rows) for three 
different charge states of the nanopore. Units: kJ/mol 







bulk water 


nanopore (-2e) 


nanopore (-5e) 


nanopore (-8e) 
















ex 






ex 


f^efe 


l^'hyd 




8.9 


-420.1 


-411.2 


1.9 


-351.3 


-349.4 


1.8 


-512.7 


-510.9 


2.0 


-696.9 


-694.9 


K+ 


22.3 


-306.9 


-284.5 


28.2 


-299.3 


-271.1 


29.7 


-454.6 


-424.9 


32.1 


-630.3 


-598.3 


Cs+ 


29.0 


-259.7 


-230.7 


58.4 


-280.1 


-221.7 


61.5 


-398.6 


-337.1 


65.2 


-576.9 


-511.7 


Na+ 


in nanopore in vacuum — > 


-4.8 


-442.1 


-446.9 


-4.8 


-1017.4 


-1022.2 


-4.8 


-1592.1 


-1596.9 


K+ 


in nanopore in vacuum — > 


15.1 


-441.3 


-426.2 


15.1 


-1014.3 


-999.2 


15.1 


-1588.6 


-1573.5 


Cs+ 


in nanopore in vacuum 


43.0 


-440.9 


-397.9 


43.0 


-1014.4 


-971.4 


43.0 


-1587.3 


-1544.3 
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TABLE ni: Free energies of water-to-hydrated-nanopore transfer of cations for three charge states of the 
nanopore. Contributions from LJ and electrostatic interactions are also Usted. Units: kJ/mol 





nanopore (-2e) 


nanopore (-5e) 


nanopore (-8e) 






A^r,d 


A^f/ 




A/^r,. 


A^r; 


A^^fe 




Na+ 


-7.0 




68.8 


61.8 


-7.1 




-92.6 


-99.7 


-6.9 




-276.8 


-283.7 


K+ 


5.9 




7.6 


13.5 


7.4 




-147.7 


-140.3 


9.8 




-323.4 


-313.6 


Cs+ 


29.4 




-20.4 


9.0 


32.5 




-138.9 


-106.4 


36.2 




-317.2 


-281.0 
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0246824682468 

simulation time (ns) 



FIG. 1: Kinetics of ion partitioning from nonequilirbium MD simulations. Number of cations partitioned 
into the negatively charged nanopore in water as a function of time in selected simulation runs. Nanopore 
carries a charge of —ne, n = 2, . . . 8 from bottom to top. Tj is the incremental time for the i-th cation to 
partition into the pore. Data for iVa+ (left), K'^ (middle), and Cs+ (right column) are shown. 
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FIG. 2: Water structure in neutral and charged nanopores. (a) Snapshots of water molecules, (b) local 
densities, and (c) orientations of water molecules in the interior of neutral and charged nanopores carrying 
a total charge of — 2e and — 5e, respectively. 6 is the angle between dipole vector and the Z axis as shown 
schematically in panel (c). 
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FIG. 3: Average orientational preferences of water molecules in the nanopore interior in the presence of a 
neutral LJ solute (K^) or a cation {K~^) for different charge states of the nanopore. Panels (b), (c), and (d) 
show snapshots from simulations. 
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FIG. 4: Free energy, Afif^^^, of water-to-hydrated-nanopore transfer of cations. LJ and electrostatic con- 
tributions to the free energy are also listed. The threshold charge density at which the transfer free energy 
becomes favorable is indicated by arrows. 
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2 3 4 5 6 7 8 

r(A) 

FIG. 5: Cation- water oxygen radial distribution functions obtained from MD simulations of cations in bulk 
water. The height of the first peak decreases significantly, indicating weaker binding of water to the ion 
with increasing cation size. Inset: Potential of mean force, or the free energy of spontaneous fluctuations 
of hydration shell waters in kT units, F{Nhy^) /kT = —\-a{p{Nhyd)], where p{Nhy^) is the probability 
distribution of number of hydration shell water molecules. 



22 



